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ABSTRACT 

We address the question how a correspondence between the particle like excitations in 
the one dimensional Hubbard model (i.e. "holons" and "spinons") and the free fermionic 
picture can be estabilished in the limit of vanishing interaction by studying the finite size 
spectrum in the framework of the Bethe Ansatz. Special attention has to be paid to the 
case of a vanishing magnetic field where the two bands of excitations in either description 
are degenerate. The interaction lifts this degeneracy. 
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The Hubbard model is one of the most studied models for interacting electrons on a one-dimensional 
lattice. Following Lieb and Wu's [Jl|] Bethe Ansatz (BA) solution many exact results have been ob- 
tained. These provide detailed understanding of the thermodynamics [0], excitation spectrum ||3j]-[§], 
finite size corrections [7|] and asymptotics of correlation functions [§, ^ which are believed to show 
the generic behaviour of systems of interacting electrons in one spatial dimension. 

The Hamiltonian is given in terms of standard fermionic creation resp. annihilation operators ^ a and 
\&j )CT of electrons with spin a at site j and the corresponding occupation numbers n 3 - )(T = ^ i0 . No- 



where u is the strength of the on-site Coulomb repulsion, fi the chemical potential and h an additional 
external magnetic field. In the following we consider the repulsive case u > with h > . 

For vanishing coupling constant u = the model simply describes two independent generations of 
free fermions. This fact allows, of course, to extract the physical properties of the system in a much 
simpler way than through the BA. Nevertheless, there have been studies of the BA solution in the 
u — > limit, mainly motivated by the desire to have a reliable test for perturbative schemes expanding 
around the free fermionic limit. These studies have mainly concentrated on the w-dependence of the 
ground state energy for the half-filled band, finding that the ground state energy can be expanded in 
an asymptotic series in u which is reproduced correctly by standard perturbation theory 

In this letter we extend the study of the u — > limit to include the behaviour of the low-lying exci- 
tations. In particular, we address the question how a correspondence between the particle like excita- 
tions of the interacting system (i.e. "holons" and "spinons") as obtained from the BA and those present 
in the free fermion spectrum can be established in this limit. It turns out that in a finite magnetic field 
there is a one to one correspondence between the excitations in either description. In absence of a 
magnetic field the single particle energies of rcomnteracting spin-| and spin-! electrons are identical. 
This is reproduced by the BA solutions where the Fermi velocities of spin- and charge excitations are 
degenerate in the limit u — > . As a consequence both pictures lead to a correct description of the 
spectrum of (at least low-lying) excitations, charge an spin degrees of freedom can be assigned to the 
two degenerate bands of excitations in an almost arbitrary way. For finite values of u this degeneracy 
is lifted. We investigate how the difference between the two velocities develops as a function of the 
interaction strength u for zero magnetic field. 

We recall that the density of charge and spin waves in the thermodynamic limit is given in terms of 
an inhomogeneous integral equation 
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where p and p^> are column vectors with entries 



and K is a 2 x 2 matrix whose elements are integral operators, namely 
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The renormalized energies of the corresponding excitations read 

e = e<® + K T * e 
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The integral operator matrix K is the transpose of K, namely 
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In these equations the kernels K\ and are given by 
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For a comparison to the free fermionic description the finite size corrections of the BA spectrum as 
calculated in [0] are particularly useful. The energies and momenta of the low lying excitations are 
given by 
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Here V denotes the diagonal matrix V = diag(t> c , v s ) of the Fermi velocities of charge and spin 
waves 
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The matrix 



'Z'cc ^cs 1 | £cc(^o) £sc(&o) 
i ^sc Z ss J \ ^ cs (Ao) £ss(Aq) 
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is given in terms of the dressed charge matrix £ which is defined by the integral equation 



£ = I + K *£ 



(13) 



where i" is the 2x2 unit matrix. The vectors 
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and the positive integers and Nf 1 characterize the excited state. Here AiV c = AA^ + AN± and 
AN S = AiVj are related to the change in particle numbers with respect to their ground state values 
thus determining charge and spin of the excitated state, respectively. D c = and D s = D± — 
are given by the number of particles moved from the left to right Fermi points at ±k F a = ±.-Kn a (n a 
are the total densities of electrons with spin a). Their values are integers or half integers subject to the 
conditions D c = (AN C + AN S ) /2 and D s = AN c /2 modulo 1. The values of N^ s are the quantum 
numbers of particle-hole excitations at the right, resp. left Fermi points. 

For vanishing u the kernels ([5]) become 5-functions and the solution of Eqs. (Q), (|5J), (jT3|) is trivial. 
However, to determine the Fermi velocities < JTT| ) and the matrix Z ( |T2| ) these solutions have to be 
taken at the boundaries k and Ao . For sin k < \ the solutions are discontinuous at these points 
and the limit u — > has to be performed after solving the integral equations. 

To see whether this situation can arise we restrict ourselves to Ao < sin k first. In this case the 
discontinuities are moved away from the boundaries entering ([TTJ) and ( |12"1 ) and we find 



(Alternatively one can express the density p s as a function of quasimomenta k = arcsin A rather than 
the rapidities A themselves, ( ]16| ) simplifies to p s (k) = 9(k — \k\)/(2ir).) From these equations we 
obtain for the total densities of the charge and spin excitations corresponding to this state 
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which allows for the identification of k and Ao in terms of the Fermi momenta through k = fcjr^ 
and arcsinAo = kpi ■ Thus we find that the condition A < sinA; is satisfied for any h > 0. The 
case of a vanishing magnetic field has to treated separately. 

The dressed energies are given by 

2u — 4cos/c if < \k\ < arcsin An 
e c (k) = { ^ ~ 1 1 " ° , (18) 

jj, — ! — 2cosfc if arcsinAo < | k\ 

, u + |-2cos£; if < \k\ < k 
e s (k) = { * 2 " . (19) 

h if k < \k\ 

From Eq. (JTTJ) we find v c = 2sm(k F ^) and v s = 2sin(/c^J . Similarly the result for the dressed 
charge matrix gives 

<i i 

1 

Now comparing the finite size corrections for the excited states (0) with this expression for the matrix 
Z and the corresponding free fermion result the two are found to agree. 
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The case of a vanishing magnetic field h — needs a special treatment. In this case we have A = oo 
and the dressed charge matrix can be expressed in terms of a single quantitity £ [0] 
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satisfying the integral equation 
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with the kernel 
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and z = sin k/u . For large z = sin k /u the quantity £(zq) entering (|T^) can be obtained using a 
perturbative scheme based on the Wiener-Hopf method JT3| ] . The result to order 1/zq reads [Q] 

£(z ) = V2(l--}-) . (24) 
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In the limit u->0we find the following dressed charge matrix 
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One might have expected that the result (EO) for Z holds even for a vanishing magnetic field as there 



is no dependence on h. This is indeed the case. "Holons" and "spinons" are certain combinations of 
spin-| and spin-| electrons. For vanishing magnetic field these combinations become arbitrary since 
spin-| and spin-| electrons have equal energies. The Fermi velocities are equal, v c = v s = 2 sin kp, 
and thus the matrix V is proportional to the unit matrix, V = 2 sin(&F) . Of physical relevance 
are only the combinations {Z~ X ) T V Z^ 1 and ZV Z T which enter expression (Q) for the excited 
states. For both choices of Z the results coincide. 

The degeneracy of the Fermi velocities of charge and spin wave exciations is lifted by the interaction. 
Using ( JTT| ) the Fermi velocities for h = can be expressed as 
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Here f(z) and g(z) are are the density p c and the derivative of the dressed energy e' c as a function of the 
variable z given in terms of the following integral equations (remember that we have uz = sin k < 1) 

1 f z " 
f(z) = + / K(z-z')f(z')dz' , 

g{z) = + / K(z- z')g(z')dz' (27) 

with the kernel K given by Eq. (fZ3J). Again, the quantities necessary to compute the Fermi velocities 
( |26[ ) for small u, i.e. z w oo , can be obtained from these equations using the Wiener-Hopf method. A 
complication is given by the explicit u-dependence of the driving terms. However, for small densities 
(i.e. uzo = sin k 1) they can be expanded up to linear order in uz. For / this results in Eq. ( |2~2"1 ) 
for the dressed charge £ (up to a factor of l/2n). In the equation for g the driving term is replaced by 
2uz. For < u sin(A; ) < 1 we finally get the following result for the Fermi velocities 
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The leading term 2uz is simply the free fermion result 2 sin k . The logarithmic corrections oc 
ulnu are probably just a consequence of the expression of the velocities in terms of z rather than 
the electron density n c . To prove this analytically the Wiener-Hopf scheme mentioned above has 
to be performed to order Zq 2 which raises questions in its quality. However, numerical solution 
of the integral equations ( |2"7| ) suggests the absence of logarithmic corrections in v a (n c , u). In Fig. 
[1] we present the Fermi velocities for a fixed value of sinfc = 0.1 which are computed from the 



numerical solution of the integral equations in comparison with Eqs. fl28P. Because of the various 



approximations which were necessary to derive Eqs. ( |2~8| ) we expect the results only do be correct up 
to the order of u. 

An interesting observation is that, in leading order, the gap between charge and spin wave excitations 
is a linear function of the interaction u 

v c -v s = — . (29) 

7T 

Fig. |2| shows the difference of the Fermi velocities as a function of the total density of particles for 
various values of u. As for Fig. |T] the data were computed from numerical solutions of the integral 
equations. As long as p, i.e. sinfc , is not to small, we find exactly the behaviour as predicted by 
Eq. (]2~9"1). For p — ► at fixed u one has z — ► which allows to solve Eqs. ([27]) by iteration. 

In this letter we have extended previous studies of the ground state properties of the one dimensional 
Hubbard model for small interaction u to the low-lying excitations. Apart from providing the possi- 
bility for a check of perturbative methods our results emphasize the importance of the "spinon-holon" 
picture for strongly correlated electrons in particular in the case of a vanishing magnetic field. 

This work has been supported in part by the Deutsche Forschungsgemeinschaft under Grant No. 
Fr 737/2-1. 
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Figures 
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Figure 1: Fermi velocities for h — and sin k Q = 0.1 as a function of u. Solid lines correspond to 
numerical solutions of the integral equations (|27|), dashed lines to the asymptotic expressions (|Z%|). 
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Figure 2: Difference of Fermi velocities for h = as a function of the total density computed from 
numerical solutions of Eqs. ([FT]). 
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